Distribution of the superconducting gap in an YNi 2 B 2 C film studied by point contact 

spectroscopy. 
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The differential resistances Rd = dV/dI(V) of point contacts between a normal metal and a c 
axis oriented YM2B2C film (T c = 15.2 K) in the superconducting (SC) state have been investigated. 
Rd{V) contains clear "gap" features connected with processes of Andreev reflection at the boundary 
between normal metal and superconductor that allow the determination of the SC gap A and its 
temperature and magnetic field dependence. A distribution of A from A m i n » 1.5 me V to A max ~ 
2.4 meV is revealed; however the critical temperature T c in all cases corresponded to that of the film. 
The value 2A max /fcBTc ~ 3.66 is close to the BCS value of 3.52, and the temperature dependence 
A(T) is BCS-like, irrespective of the actual A value. It is supposed that the distribution of A can 
be attributed to a gap anisotropy or to a multiband nature of the SC state in YM2B2C, rather than 
to the presence of nodes in the gap. 

PACS numbers: 74.70.Dd, 74.45. +c, 74.50.+r 



INTRODUCTION 



During the last ten years, the class of nickel borocar- 
bide compounds i?Ni 2 B 2 C, where R = Sc, Y, Lu, Dy, Ho, 
Er, Tm and other rare-earth ions, has attracted consid- 
erable attention in the superconducting community 0. 
This interest is caused by the unusual properties of these 
materials related to the rich interplay of superconduc- 
tivity and magnetism, when R is a magnetic rare-earth 
ion. In this context, it is remarkable that the supercon- 
ducting transition temperature T c can be higher than the 
Neel temperature Tn as for i?=Tm and Er, but also com- 
parable with Tn as for i?=Ho and even lower than Tn as 
for R=Dy. Non-magnetic compounds with R—Sc, Y and 
Lu have the highest T c values of 15 K, 15.5 K and 16.5 K, 
respectively. They are type II superconductors and have 
three-dimensional, practically isotropic electronic struc- 
ture, while their crystal structure comprises alternating 
RC and N2B2 layers. Noteworthy recent experiments in- 
cluding tunnelling 0] and point contact measurements 
0, 0] on bulk single crystals have shown that the un- 
usual superconducting properties of nickel borocarbides 
exhibit strong anisotropy. In particular, a hybrid (s + g) 
wave superconducting state has been proposed 0, |(J for 
the compounds with R = Y and Lu, characterized by 
four point nodes of the superconducting order parameter 
in the [±1,0,0] and [0,±1,0] directions. This conclusion 
was based on investigations of the angular dependence of 
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the thermal conductivity and specific heat || under 
an applied external magnetic field of single crystals at 
low temperatures. Point-contact measurements Q also 
revealed a strong gap anisotropy. In the [100] direction, 
not only the gap is 4.5 times smaller, but also the transi- 
tion temperature was found to be half (!) that in the [001] 
direction. It is also important to note that a thorough 
analysis of upper critical field H C 2(T) data for YM2B2C 
and LUM2B2C single crystals clearly showed that super- 
conductivity in these compounds is strongly affected by 
their multiband nature |9(- In order to gain further in- 
sight into the interplay of the symmetry of superconduct- 
ing order parameters, their anisotropy and multiband 
nature, additional studies are necessary. In the present 
work, direct measurements of the gap amplitude distri- 
bution have been performed on epitaxial c axis oriented 
YNi 2 B 2 C thin film using point-contact spectroscopy [Io| . 

EXPERIMENTAL TECHNIQUE 

In this paper, we present the results of point contact 
measurements made on a high quality YM2B2C film [ill ]. 
The film was prepared by pulsed laser deposition on an 
MgO single crystal substrate. The SC transition temper- 
ature T c = 15.2 K (with transition width AT C = 0.3 K) 
was determined from resistive measurements. The resid- 
ual resistance ratio R(30QK) / R(T > T c ) of the film was 
about 8. X-ray data indicate that the film consists of 
small crystallites having their primary orientation with 
the c axis perpendicular to the substrate surface. 

The measurement cell was placed in a flow cryostat, 
enabling measurement at a series of temperatures up to 
T c and higher. Point contacts were created by touching 
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FIG. 1: (a) Differential resistance Rd(V) of point contact 
YNi2B2C-Ag with normal state resistance Rn = 3.5J1, mea- 
sured at different temperatures, (b) Normalized Rd (solid) 
curves for this contact and those values calculated by BTK 
theory (dashed). The theoretical curve at 4.2 K is calculated 
with A « 1.67 meV, Z » 0.43 and F « 0.55 meV. The tem- 
perature dependence A(T) is shown in Fig.[3Jb). Curves are 
offset vertically from the lower curve for clarity. 



FIG. 2: (a) Differential resistance Rd(V) of point contact 
YN12B2C-CU with normal state resistance Rn ~ 4.5f2, mea- 
sured at different temperatures, (b) Normalized Rd (solid) 
curves for this contact and those values calculated by BTK 
theory (dashed). The theoretical curve at 4.5 K is calculated 
with A « 2.24 meV, Z w 0.5 and F tti 0.1 meV. The tempera- 
ture dependence A(T) is shown in Fig.^b). Curves are offset 
vertically from the lower curve for clarity. 



sharpened Cu and Ag wires to the film surface directly 
within the cryostat at liquid helium temperature. An 
unfortunate disadvantage of this 'needle-anvil' method 
is the sensitivity of the contacts to mechanical vibrations 
and changes in temperature. As a result, only around one 
quarter of the prepared contacts withstood measurement 
above liquid helium temperature 

Using a technique of synchronous detection of weak 
(10-30 /iV) modulating signal harmonics, the first har- 
monic of the modulating signal V\ [proportional to the 
differential resistance R<i(V) — dV / dl (V)] was recorded 
as a function of the bias voltage V . 



EXPERIMENTAL RESULTS 

A series of Rd{V) measurements with clear double min- 
imum features near V — connected with the Andreev 
reflection of electrons at the superconductor - normal 



metal (S-N) interface are shown in Figs.^a) and^a). 

For voltages exceeding the gap value, all curves exhibit 
maxima caused, most likely, by the suppression of super- 
conductivity in (part of) the contact. Such behavior tes- 
tifies to a deviation from the spectral (ballistic) regime 
(i.e., the inelastic electron mean free path becomes less 
than or about the size of the contact) with increasing 
bias voltage at the contact (see 0, Chapter 3). 

The size of the gap was calculated using the modi- 
fied BTK theory (see |lfj, Chapter 3.7). The presence of 
maxima in Rd(V) prevented the coincidence of the calcu- 
lated and measured curves throughout the voltage range; 
therefore an adjustment was carried out so that both 
curves coincided in the region of the gap (—A, +A) and 
at the higher voltages beyond the maxima. Preliminary, 
unadjusted curves were normalized to Rd{V, T > T c ). 
During calculation, the parameters Z and T were fixed 
or varied in small limits up to 5-10%. The results are 
shown in Figs.^b) andj^b). We note also that the con- 
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FIG. 3: (a) The temperature dependence of the gap A(T) 
obtained from the BTK fit for five different contacts (open 
symbols). The theoretical BCS curve is drawn to coincide 
with the average value of the five contacts (closed symbols), 
(b) The temperature dependence A(T) normalized to A(T = 
4.5 K) for minimal (open squares), maximal (open circles) and 
average (closed circles) gap values, along with the theoretical 
BCS curve. 




FIG. 4: Distribution of A values obtained from measure- 
ment of 29 contacts with different counter electrodes Ag (15 
contacts, dark bars) and Cu (14 contacts, gray bars). The 
arrow specifies the value corresponding to the BCS ratio 
2A /fc B Tc = 3.52, with T c = 15.2 K. 




tacts shown in Figs.nand[21had, respectively, the mini- 
mal and maximal gap values among all the contacts that 
withstood low temperature measurement. 

In Fig.Efa), the dependences A(T) are shown for a 
number of different contacts, for each of which a series 
of Rd{V) curves at different temperatures was measured. 
For all these contacts, the SC minimum in Rd(V) disap- 
peared at a temperature interval of 14-16 K, correspond- 
ing to the T c value (15.2 K) of the film. 

In Fig.[3{b), the temperature dependences for the max- 
imal and minimal A(T = 4.5 K) are shown, normalized to 
these values (2.24 meV and 1.67 meV, respectively). Also 
the normalized A(T) temperature dependence obtained 
by averaging all five contacts of Fig.[3Ja), and the theo- 
retical BCS dependence are depicted. We reiterate that 
the complete data (measured, normalized and calculated 
spectra) for the smallest and largest gap are presented in 
Figs.HlandEl 

In total, 29 curves, on which the structure with two 



FIG. 5: Distribution of contact resistances Ro = Rd(V = 0, 
T — 4.2 K) for the same contacts as in Fig.^] Inset: gap 
values of these contacts. Solid circles designate those contacts 
for which a series of temperature measurements were carried 
out, with results shown on Fig.^a). 

minima is clear enough to allow an unequivocal determi- 
nation of A using the standard procedure of adjustment 
of theoretical curves described above, have been selected 
from numerous measured spectra. From the results of the 
calculation, the histogram in Fig.^was constructed. Ir- 
respective of the counter electrode used, A is distributed 
in an interval between 1.5 meV and 2.4 meV. 

Half of the selected contacts had contact resistances 
between 2 and 3 f2, as seen on the distribution of N(R) 
values in Fig.[5j The inset shows the A values obtained 
for contacts of different normal state resistance. Solid 
circles designate those contacts for which a series of tern- 
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FIG. 6: (a) Differential resistance Rd(V) of the same point 
contact YN12B2C-CU measured for different Ro — Rd(V = 0, 
T — 4.2 K) values as a result of consecutive contact short- 
circuits by a current pulse, (b), (c), (d) The variation in the 
values A, F and Z with Ro. 



perature measurements were carried out and supercon- 
ductivity confirmed to disappear between 15 K and 16 K, 
corresponding to the T c value of the film, independent of 
the A value. 

In Fig.Efa), Rd(V) is shown for contacts produced by 
consecutive electric short-circuit of the same contact from 
Rn = 20 fi down to Rn = 4 fl. It is seen that as the con- 
tact resistance is reduced, humps start to appear to a 
greater degree, and the intensity of the local maximum 
at zero bias increases. The latter is possible to explain, 
proceeding from ^3] where it is shown that in the dif- 
fusive regime Rd(V) looks like the presence of a barrier 
with Z= 0.55 in the ballistic regime. In other words, 
for a clean S-N boundary without a barrier (Z = 0) in 
the diffusive regime, it is possible to observe a maximum 
in Rd(V) at zero bias, similar to the ballistic case with 
Z = 0.55. In reality, from Fig.[f)fc) it is seen that the 
parameter Z grows with the reduction of contact resis- 
tance, coming nearer to the value 0.55. It is necessary 
to note that the A value for these contacts is practi- 
cally independent of their resistance, remaining around 
2.2 ± 0.1 meV. On the other hand, the parameter T grows 
strongly with resistance, as observed experimentally for 
many superconductors of different sorts (see Fig. 12.14 in 
jlOj). although the reason for this remains under discus- 
sion. 



DISCUSSION 

The presented point contact measurements show a dis- 
tribution of the SC gap in the electronic density of states 



of an YNI2B2C film. Important and essential is the fact 
that irrespective of the size of the gap, ranging from 
1.5 meV to 2.4 meV, the observed critical temperature 
corresponds to the T c of the film. This excludes the pos- 
sibility of a distribution of the gap on the film surface ow- 
ing to its non-uniform structure. Another possible source 
of a gap variation is the proximity effect resulting in gap 
suppression at the surface. However, we assume that in 
this case the gap value should correlate inversely with 
the point contact resistance, in other words with the size 
of the contact: the smaller the size, the lower the gap 
value measured. However, this is not the case. The gap 
value varies even for contacts with nearly equal resistance 
[Fig-EImset)], and has a tendency, if any, to rise with in- 
creasing resistance [Fig.EJb)]. 

Hence, it is natural to first assume the presence of a 
gap anisotropy. Unfortunately, the use of a film does not 
allow measurements to be made along different crystal- 
lographic directions. However, as mentioned previously, 
the film is predominantly c axis oriented, and consists of 
small crystallites. Thus it is reasonable to expect that 
with contacts prepared by touching a film plane with a 
needle, we will mainly register a gap along the c axis. 
Yet, taking into account that the contribution to point 
contact conductivity result from some finite solid angle, 
the point contact axis is not well fixed along the c axis 
and with the possible misorientation of some crystallites 
our data can be interpreted as the observation of the gap 
anisotropy with A max /A m ; n w 1.6. We also note that, 
as seen from Fig.0J the ratio 2A max /fceT c is close to the 
standard BCS value of 3.52. 

As to the presence of a hybrid gap with s + g wave 
order parameter symmetry and point nodes in certain 
directions, we never observed gaps less than 1.5 meV. 
The results of point contact measurements in giv- 
ing A 0.45 meV along the a axis, give rise to some 
questions. Firstly, the dependence A(T) for the a direc- 
tions looks at least atypical; secondly, the critical tem- 
perature in this direction was found to be two times less 
than in the bulk; and thirdly, even for the maximal gap, 
2A max //cBT c ss 2.86, much lower than the BCS value. 
This raises the doubt of whether to take into account 
that nickel borocarbides are more likely superconductors 
with strong coupling. We note that the single crystals 
investigated in |3j had a significantly lower T c w 14.6 K, 
in comparison with the value 15.5 K usually quoted in 
the literature. Thus, an insufficient sample quality is ap- 
parently responsible for the reported " anisotropy" of the 
SC order parameter in ||. 

Considering the role of disorder in anisotropic super- 
conductors, nonmagnetic impurity scattering has pro- 
found effects on hybrid s + g wave superconductors with 
point nodes. It is mentioned in |6| that a disorder- 
induced quasiparticle energy gap opens up even for in- 
finitesimal scattering rates. This may be a reason why 
we did not observe any nodal features in the SC gap. 
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However, impurity scattering also decreases the critical 
temperature, which is not observed in our investigations. 
Thus, according to 0, a disorder- induced energy gap of 
about 10 K w 1 meV corresponds to a T c reduction to 
13K. 

Concerning the effect of the pressure produced by 
point-contact formation on the gap distribution, we ex- 
pect that under such a pressure the critical temperature 
should also vary in the same way as the gap. However, 
no such variation is observed, suggesting that there is no 
measurable effect due to pressure in this experiment. 

Let's consider the question of the realization of multi- 
band superconductivity in i?Ni 2 B 2 C (R = Y, Lu), con- 
sidered in 0. We cannot distinguish whether the ob- 
served gap distribution is a consequence of anisotropy or 
a multiband electronic structure. In this connection, we 
note that the authors of Ref. jj] , investigating point con- 
tact spectra of a LuNi 2 B 2 C single crystal, found that the 
calculated theoretical curves better describe the experi- 
mental data when a two-gap adjustment is used, although 
in the Rd{V) curves, only one gap minimum is visible. In 
so doing, the relation between the maximal and minimal 
gap was found to be around 1.75 in the c direction and 
1.5 in the ab plane. This is close to the value found by 
us, A max /A min rs 1.6. To this it is possible to add that 
our preliminary data after analysis of about 60 Rd(V) 
spectra measured on a cleaved YNi 2 B 2 C single crystal 
also show a gap distribution. In this case for 80% of the 
contacts the gap is distributed in the same range from 
1.5 meV to 2.4 meV, while for 20% of the contacts the 
gap ranged more from 2.5 meV to 3.0 meV, and one can 
discern maxima in the entire gap distribution at 2.0 meV 
and 2.4 meV. Interestingly, to fit the experimental dV/dl 
curve for the c direction in LuNi 2 B 2 C in Ref. |4j (Fig. 13), 
the authors used a gap distribution with two peaks also 
around 2.0 meV and 2.5 meV. 

A prominent feature of a multiband superconductor 
is also the specific behavior of the excess current in the 
point contacts. In |l3| it was revealed that in the two- 
band superconductor MgB 2 , the excess current has a 
strongly pronounced positive curvature in a magnetic 
field. In Fig.Q the behavior of the SC gap and the ex- 
cess current in a magnetic field for one of the YNi 2 B 2 C 
contacts is shown. It can be seen that, while the be- 
havior of the gap is similar to that of the pair potential 
of a type II superconductor in the vortex state and has 
negative curvature [Fig.0a), closed squares and dotted 
curve] , the excess current has a strong positive curvature 
[Fig-Efb)]. A similar behavior of the gap and the excess 
current in a magnetic field was observed by us for several 
of the investigated contacts (Rd(V) have been processed 
for four contacts). Another consideration is that a small 
gap has an almost linear decrease under a magnetic field 
[Fig-Efa), open squares]. This is similar to the theoreti- 
cal behavior of the small gap in MgB 2 calculated in . 
Thus, our observation of the excess current and small gap 
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FIG. 7: (a) Calculated magnetic field dependences of A 
(closed squares) and V (closed circles) at Z — 0.58 for an 
YNi 2 B 2 C-Cu point contact (R = 2.8 £1, T = 4.2 K). Open 
squares show the magnetic field behavior of A for another 
contact with small gap value. The dotted curve shows the 
behavior oc (l-B/B c2 ) 1/2 of the pair potential of a type II su- 
perconductor in the vortex state according to Abrikosov's the- 
ory. Dashed lines connect circles and open squares as guides 
to the eye. Inset: reduced experimental Rd(V) in zero field 
and in magnetic fields of 1 T and 5 T. (b) Behavior of the ex- 
cess current in a magnetic field. The inset shows symmetrized 
experimental and theoretical R^ (V) in zero field. The ex- 
cess current was determined as the integral of the calculated 
iZj^V) after subtraction of the normal state conductivity 
(horizontal dashed line). Error bars were determined as the 
integral of the difference between theoretical and experimen- 
tal R^ (V) curves. The greatest difference arises due to the 
presence of minima in the experimental RJ 1 (V) curves at bi- 
ases exceeding the gap energy. 

behavior under a magnetic field may be taken in support 
of a multiband SC state in YNi 2 B 2 C. 

CONCLUSION 

Studies of the SC gap A and its temperature depen- 
dence A(T) in YNi 2 B 2 C have shown a distribution of 
A values from A m i n w 1.5 meV up to A max w 2.4 meV, 
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however in all cases with a BCS-like A(T) dependence. 
The observed 2A max /fcBT c rj 3.66 is also close to the 
BCS value of 3.52, and the critical temperature in all 
cases corresponded to that of the film, thus excluding 
the possibility that the observations result from hetero- 
geneity in the properties of the film. It is most probable 
that the distribution of the gap can be attributed to gap 
anisotropy, or to the multiband nature of the SC state in 
YNi2B2C. The positive curvature in the behavior of the 
excess current under applied magnetic field (Fig.Efb)) 
tends to support the multiband scenario. Proceeding 
from the stated results, and also preliminary investiga- 
tions of YNi2B2C single crystals, the presence of point 
nodes in the SC order parameter of YM2B2C is consid- 
ered to be unlikely. The analogous conclusion is reached 
also in a recent work on the basis of the observation 
of a sudden disappearance of de Haas-van Alphen oscil- 
lations at the transition to the mixed state. Finally, in a 
recent preprint |l6j . the authors of Ref.Q have modified 
their earlier claims and come to the conclusion from the 
analysis of the field variation of their dV/ dl spectra that 
the unconventional gap function observed in YNi2B2C 
could originate from multiband superconductivity. 
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